Search results for "numerical [Methods]"
showing 10 items of 500 documents
Random wave run-up with a physically-based Lagrangian shoreline model
2014
Abstract In the present paper the run-up of random waves was calculated by means of a numerical method. In situ measurements based on a video imaging technique have been used for the validation of the present numerical model. The on-site run-up measurements have been carried out at Lido Signorino beach, near Marsala, Italy,along a transect, normal to the shore. A video camera and a linear array of rods have been used to obtain field data. Numerical simulations with a 1DH Boussinesq-type of model for breaking waves which takes into account the wave run-up by means of a Lagrangian shoreline model have been carried out. In such simulations random waves of given spectrum have been propagated in…
Testing steady-state analysis of single-ring and square pressure infiltrometer data
2016
Testing reliability of the saturated soil hydraulic conductivity, Ks, estimated by applying the steady-state single-ring (SR) model to the quasi steady-state infiltration rates obtained with a single-ring pressure infiltrometer (PI) increases confidence in the estimated Ks values. Determining a means to estimate Ks from infiltration data collected with a square infiltrometer allows the use of sources of different shapes. Using numerically simulated infiltration rates for six homogeneous soils ranging in texture from sand to silty clay loam, this investigation suggested an overall good performance of the SR model, with estimated Ks values differing by not more than 25% from the true values f…
Saturation excess runoff numerical simulation
2005
Saturation excess runoff is a relevant process which needs additional experimental and modeling efforts. This work is focused on its numerical modeling. The final objective is the successive interpretation of ongoing experimental monitoring results in two watersheds in different areas of Italy where the saturation excess runoff formation mechanism seems to be important. The numerical solution of the two-dimensional Richards’ equation allows the evaluation of the sensitivity to the various influent parameters : rainfall intensity, soil properties, depth and initial water content, slope and hillslope length. Also the subsurface flow is simulated at the same time, allowing the evaluation of th…
4D paleoenvironmental evolution of the Early Triassic Sonoma Foreland Basin (western USA)
2017
In the wake of the Mesozoic, the Early Triassic (~251.95 Ma) corresponds to the aftermath of the most severe mass extinction of the Phanerozoic: the end-Permian crisis, when life was nearly obliterated (e.g., 90% of marine species disappeared). Consequences of this mass extinction are thought to have prevailed for several millions of years, implying a delayed recovery lasting the whole Early Triassic, if not more. Several paradigms have been established and associated to a delayed biotic recovery scenario expected to have resulted from harsh and deleterious paleoenvironments. These paradigms include a global anoxia in the marine realm, a “Lilliput” effect, and the presence of “disaster” tax…
Bill2d - a software package for classical two-dimensional Hamiltonian systems
2015
Abstract We present Bill2d , a modern and efficient C++ package for classical simulations of two-dimensional Hamiltonian systems. Bill2d can be used for various billiard and diffusion problems with one or more charged particles with interactions, different external potentials, an external magnetic field, periodic and open boundaries, etc. The software package can also calculate many key quantities in complex systems such as Poincare sections, survival probabilities, and diffusion coefficients. While aiming at a large class of applicable systems, the code also strives for ease-of-use, efficiency, and modularity for the implementation of additional features. The package comes along with a use…
Analytical evaluation of structural response for stationary multicorrelated input
1990
Abstract An analytical procedure is presented which can drastically reduce computational effort in the evaluation of the spectral moments of an elastic linear multi-degree-of-freedom system subjected to a stationary multicorrelated input process. The reduction in computer time is possible since the cross-spectral moments of two oscillators can be obtained in recursive manner as a linear combination of the spectral moment of each oscillator taken separately, which is evaluated by means of a very fast numerical technique.
Some efficient algorithms for the solution of a single nonlinear equation
1981
High order methods for the numerical solution of nonlinear scalar equations are proposed which are more efficient than known procedures, and a unified approach to various methods suggested in literature is given.
Multi-Phase epidemic model by a Markov chain
2008
Abstract In this paper we propose a continuous-time Markov chain to describe the spread of an infective and non-mortal disease into a community numerically limited and subjected to an external infection. We make a numerical simulation that shows tendencies for recurring epidemic outbreaks and for fade-out or extinction of the infection.
Simulation of BSDEs with jumps by Wiener Chaos Expansion
2016
International audience; We present an algorithm to solve BSDEs with jumps based on Wiener Chaos Expansion and Picard's iterations. This paper extends the results given in Briand-Labart (2014) to the case of BSDEs with jumps. We get a forward scheme where the conditional expectations are easily computed thanks to chaos decomposition formulas. Concerning the error, we derive explicit bounds with respect to the number of chaos, the discretization time step and the number of Monte Carlo simulations. We also present numerical experiments. We obtain very encouraging results in terms of speed and accuracy.
Newton algorithm for Hamiltonian characterization in quantum control
2014
We propose a Newton algorithm to characterize the Hamiltonian of a quantum system interacting with a given laser field. The algorithm is based on the assumption that the evolution operator of the system is perfectly known at a fixed time. The computational scheme uses the Crank-Nicholson approximation to explicitly determine the derivatives of the propagator with respect to the Hamiltonians of the system. In order to globalize this algorithm, we use a continuation method that improves its convergence properties. This technique is applied to a two-level quantum system and to a molecular one with a double-well potential. The numerical tests show that accurate estimates of the unknown paramete…